Long-range entanglement is necessary for a topological storage of quantum 

information 
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We ask whether the ground state entanglement of a quantum many-body system imposes a fun- 
damental constraint on its information storage capacity. We show that the amount of long-range 
entanglement gives a rigorous upper bound to this problem, even without invoking the properties 
of the parent Hamiltonian. In particular, we show that log TV < 27 for a large class of topologi- 
cally ordered systems on a torus, where A'^ is the number of topologically protected states and 7 
is the topological entanglement entropy. We also obtain a necessary condition for a quantum error 
correcting code to have a nonvanishing rate: its average entanglement entropy over the subsystems 
smaller than the code distance must satisfy a strict volume law. 



Entanglement entropy is a canonical measure for quan- 
tifying entanglement in a bipartite pure state. 1 There 
has been a recent surge of interest in studying entangle- 
ment entropy in a ground state of a quantum many-body 
system. [5^ One of the motivations behind these studies is 
that the entanglement entropy is a useful probe for de- 
tecting the phase of the quantum many-body system. For 
example, entanglement entropy in one dimensional crit- 
ical systems follows a universal logarithmic scaling law, 
and its prefactor is related to the conformal charge of the 
theory. [3] In two spatial dimensions, the quantum dimen- 
sion of the topological quantum field theory describing 
the low-energy physics can be inferred from a constant 
subcorrection term of the entanglement entropy. ^^ 

Another important motivation comes from the nu- 
merical simulation of quantum many-body systems. 
Classes of variational ansatz such as the matrix 
product states(MPS),[71 [5] projected entangled pair 
states(PEPS),i9J and the multi-scale entanglement renor- 
malization ansatz (MERA) [lO] have certain entropy scal- 
ing laws. Since these variational states reproduce the en- 
tanglement scaling of gapped/critical systems, they are 
suitable for efficiently simulating the ground state prop- 
erties of the quantum many-body systems. In particular, 
a rigorous argument can be made for one-dimensional 
gapped systems, where an explicit scaling relation be- 
tween the entanglement entropy and the MPS bond di- 
mension is known. pTM 13) 

Recently another possibility has been explored by sev- 
eral authors. [T3HIS] Their approaches differ from the con- 
ventional ones in several aspects. For one thing, they ex- 
pHcitly use the entanglement entropy over multiple sub- 
systems simultaneously. Moreover, the main objective of 
these works is not necessarily focused on obtaining an or- 
der parameter of the phase. Rather, they are interested 
in identifying a hidden structure of the quantum many- 
body phase that protects their universal properties. 

This hidden structure by in large can be attributed 
to the structure of states that are approximately condi- 
tionally independent. A tripartite state pabc is condi- 



tionally independent if its conditional mutual informa- 
tion I{A : C\B) = S{AB) + S{BC) ~ S{B) ~ S{ABC) is 
equal to 0, where S{A) = — Tr(p^ log pa) is the entangle- 
ment entropy of the subsystem A. Similarly, a tripartite 
state is approximately conditionally independent if its 
conditional mutual information has a nonzero value that 
is close to 0. At least for gapped quantum many-body 
systems, there is a good reason to expect such subsys- 
tems to appear quite naturally. [14] The power of these 
approaches lies on the fact that their argument only re- 
lies on the generic properties of the ground state alone. 
Moreover, the structure that arises from the vanishing 
conditional mutual information is manifestly nonlocal: it 
gives a nontrivial constraint between the density matrices 
over large regions. [TTHTn] 

In this Letter, we explore this structure further in the 
context of finding a fundamental limit on the informa- 
tion storage capacity of a quantum many-body system. 
There has been a number of results in recent years, where 
information storage tradeoff bounds for local quantum 
codes have been obtained. [20h24J A local quantum code 
refers to a set of degenerate ground states stabilized by a 
sum of geometrically local commuting projectors. Impor- 
tant examples include quantum double and Levin- Wen 
model. pSl [2S] Bravyi et al. showed that the following 
bound holds for such systems 
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where k is the number of qubits, d is the code distance, 
D is the spatial dimension of the lattice, and n is the 
number of particles. [21] 

Local quantum codes cover a rich array of systems. In- 
deed, it has been conjectured that Levin- Wen model can 
describe any non-chiral topologically ordered systems in 
two spatial dimensions. |26j However, none of these mod- 
els give rise to the chiral gapless boundary excitations. 
Hence, the tradeoff bound cannot be applied to physi- 
cal systems exhibiting the integer quantum Hall effect or 
the fractional quantum Hall effect. [27 | [28 ] Furthermore, 
there are models that are topologically ordered yet do not 



necessarily have a local commuting parent hamiltonian. 
Examples include certain quantum dimer models. |29M30j 
parent hamiltonian of PEPS,[3T1 132] and Kitaev's hon- 
eycomb model. 33] 

To deal with such systems, we approach this problem 
by focusing only on the properties of the state. More pre- 
cisely, we obtain an upper bound on the number of topo- 
logically protected qubits in terms of a linear combina- 
tion of entanglement entropies over some local regions of 
a single wavefunction. Entanglement entropies are later 
related to the code parameters. At first, the preceding 
claim might seem oxymoron to some of the readers. How 
can one bound a number of degenerate ground states 
without knowing the parent hamiltonian? Further, how 
can one even define an entanglement entropy of a subsys- 
tem if there are multiple states? Answers to both of these 
questions lie on an important difference between a classi- 
cal error correcting code and a quantum error correcting 
code. 

The key property of the quantum error correcting code 
that sets apart from its classical counterpart is the lo- 
cal indistinguishability property. While a more gen- 
eral version of this definition exists. [MH5B] it will be 
convenient to use the concept of topological quantum 
order(TQO)[37] for our purpose. Following Ref,37, we 
say that a set of states {|i/'i)}i=i,... ,Ar satisfies a TQO 
condition with (r, e)-error if 



(2) 



holds for any (j) that is restricted to a ball of radius r, 
where || • • • || is the operator norm. If the approximation 
parameters r and e are obvious from the context, we shall 
simply say that the states are locally indistinguishable. 
An important consequence of EqJ2] is that the reduced 
density matrices of the locally indistinguishable states 
are close to each other in a trace distance. Therefore, 
one can unambiguously define the entanglement entropy 
of the aforementioned set of states up to a small error, 
so long as the subsystem can be contained in a ball of 
radius r. We shall call such subsystems to be local. [59] 

Roughly speaking, each of the qubits in the quantum 
error correcting code is entangled to each other so that 
the information can be stored nonlocally. As the number 
of encoded qubits increase, more entanglement is neces- 
sary to distribute the information sufficiently nonlocally. 
Hence, there has to be an inherent limit on the number 
of protected qubits if the entanglement across different 
subsystems are bounded. 

In this Letter, we shall set this intuition on a rigorous 
ground. Our approach is purely information-theoretic: 
we only use the strong subadditivity of entropy as our 
main technical tool.[33 However, we believe it would be 
more instructive to discuss some of the related results 
to understand the motivation behind our construction. 
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Strong subadditivity of entropy asserts that the condi- 
tional mutual information I {A : C\B) of a tripartite 
quantum state Pabc is nonnegative. While this state- 
ment is true for any quantum states, a special structure 
arises when the inequality is satisfied with an equality. In 
particular, Petz[T7] showed that a tripartite state pabc 
can be reconstructed from paBtPb and pBc if the con- 
ditional mutual information vanishes: 



PABC = PabPb PbcPb Pab- 



(3) 



Recursively applying Petz's result, Hastings and 
Poulin showed that quantum systems that saturate the 
equality condition for certain subsystems can be com- 
pletely reconstructed from the local reduced density ma- 
trices alone. [H] A simple corollary of their result is that 
there cannot be another distinct state that is locally in- 
distinguishable from the original state. To understand 
why, suppose there are two distinct states pi , p2 that are 
locally indistinguishable. Since both states have the same 
local reduced density matrices, the conditional mutual 
information computed from these local density matrices 
must be both 0. Due to the same reason, a recursive ap- 
plication of Eqj3] must yield a same recovered state for 
both pi and p2- Since pi and p2 were assumed to be 
distinct, this contradicts the original assumption. While 
this argument can be applied quite generally, it is not 
robust against small deviations of the original assump- 
tions: if the conditional mutual information is only ap- 
proximately equal to 0, EqlS] cannot be used anymore. 
Also, Hastings and Poulin's result applies only to tree 
graphs. We shall discuss how these issues can be circum- 
vented by giving several concrete examples. 

-An inequality between topological entanglement en- 
tropy and topological degeneracy 

We first discuss a class of models in two spatial di- 
mensions that satisfy a strict form of area law. Loosely 
speaking, we say a model satisfies a strict form of area 
law if the nonuniversal contribution to the entanglement 
entropy can be canceled out by taking a judicious choice 
of a linear combination. The entanglement entropy of 
these models can be written as 



S{A) = a\dA\-bo{A)-f + 0{e 
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where a is a nonuniversal constant, \dA\ is the number 
of particles along the boundary of A, ^ is the correlation 
length, and 7 is the topological entanglement entropy. [H 
1^ bQ{A) is the number of connected of components of the 
boundary of A. The correlation length will be equal to 
for a fixed-point wavefunction of some renormalization 
group flow, but otherwise would remain nonzero. 

Suppose we have N states {|V'i)}i=i.--- n with 
(cL,0(e--^/«))-TQO condition that satisfy Eq|4|for a lo- 
cal subsystem. Here we have assumed our system to be 
on a L X L torus. The specific value of the numerical 
constant < c < 1 is irrelevant for the analysis, as long 



as it is sufRciently close to 1. Our first main result gives 
a rigorous lower bound of 7 in terms of N, up to a small 
correction that vanishes in the thermodynamic limit. 



log TV < 27 + C'(L^e"-^/«) 



(5) 



The idea for proving Eq[5] is to apply the Markov en- 
tropy decomposition(MED) to a maximally mixed state 
over the N states. [13] More precisely, consider a se- 
quence of subsystems Ai,Bi,Ci, i = 1 , ■ • • ,n such that 
(i) AiBiCi — Ai^iBi^i, (ii) AiBi and BiCi are locals 
and (iii) AnBnCn is the entire system. For such choice 
of subsystems, following linear combination of entangle- 
ment entropy is nonnegative due to the strong subaddi- 
tivity of entropy: 

n n 

Y, HA^ ■■ C^m = SiAiBi) + Y, S{B,a) ~ S{Bi) 

1=1 i=l 

~ ^[AnB„Cn)- 

By choosing the global state to be a uniform mix- 
ture of the A^ locally indistinguishable states, i.e. 
Xli=i W IV'i) (V^ilj we arrive at the following bound: 



logiV < SiAiBi) + Y S{B,C,) - S{B, 



(6) 



j=i 



Since AiBi, BiCi, and B^ are all local, their entangle- 
ment entropy can be replaced by an entanglement en- 
tropy of one of the states \ipj) with a small correction. 
The correction term can be estimated by using Fannes 
inequality [39] which holds for any quantum states p,a 
supported on a d-dimcnsional Hilbert space: 

\S{p)~S{a)\<e\ogd-e\oge, 6:=|p-a|i, (7) 

where | • • • 1 1 is the trace norm. EqJ5] can be derived by 
choosing an appropriate set of subsystems such that the 
boundary contributions cancel out, while the universal 
term survives. One choice of such subsystems is depicted 
in FIGdl 

While we have assumed that the corrections to the area 
law and the local indistinguishability property decay ex- 
ponentially in the system size, the same argument can 
be carried out even for an algebraically decaying func- 
tion with a large enough exponent, see the Appendix. 
Eq(5] confirms the intuition that an amount of long-range 
entanglement limits the topological ground state degen- 
eracy. For two-dimensional models supporting anyons, 
this result is already known. The ground state degener- 
acy of a topologically ordered system is rip, where Up is 
the number of particle types. [301 On the other hand, the 
topological entanglement entropy is related to the total 
quantum dimension of the system: 
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FIG. 1: Each of the diagrams represent the subsystems 
Ai,Bi,Ci, i = 1, 2, 3 with the constraint that AiBid is equal 
to Ai+iBi^i. Due to this construction, the entanglement en- 
tropy of A3B3C3 is the only nonlocal contribution to the sum 
"^i^iH^i ■ Ci\Bi). The entanglement entropy of AsBsC-j 
becomes log A'^ for a maximally mixed state over A'^ locally 
indistinguishable states. The rest of the contributions can be 
computed from the formula for the entanglement entropy of 
local subsystems, i.e. Eq[4] 



where da is the quantum dimension of a particle with a 
topological charge a. [3 [6] Since non-Abelian charges can 
have a quantum dimension that is strictly larger than 1, 
topological entanglement entropy has to be always larger 
or equal to |logiV.[S] However, we emphasize that we 
did not assume anything about the Hamiltonian at all. 
Eq[5] was derived only from the property of the states. 
Our result shows that even the topological ground state 
degeneracy - typically thought as a property inherited 
from the Hamiltonian - is already strongly constrained by 
the structure of the state alone. Admittedly our premise 
about the form of the entanglement entropy is quite re- 
strictive. However, we note that there are number of 
models whose ground state entanglement entropy can be 
either analytically computed [51 [51 [^ S5", "52", ^T or com- 
puted with a good numerical precision ,42.-.46, . 

- Higher dimensional systems 

The preceding argument is based on the observation 
that the entanglement entropy of the subsystem can be 
decomposed into two parts. The first part is the local 
entanglement, which can be canceled out by making a 
judicious choice of subsystems. More precisely, we expect 
the entanglement entropy of a subsystem A to have the 
following form: 

S{A)=Slocal{A) + Stopo{A), 

where Siocaii^) is an entanglement that can be decom- 
posed into strictly local contributions. [17] The rest of 
the contributions, including the finite-size effects, are 
included in Stopo{A). For such systems, the maximal 
number of topologically protected states depends on the 
scaling law of StopoiA). Alternatively, a quantum error- 
correcting code with an extensive ground state degen- 
eracy would imply that there is an extensive subcor- 
rection term for the entanglement entropy. Interest- 
ing examples include Chamon's model and Haah's cubic 
code, which are known to have ground state degeneracies 
that increase as A^ = 2^^^^ for certain choices of L.|481- 
150] Therefore, the entanglement entropy of these mod- 
els must have a subcorrection term that grows as il{L). 



Also, one may argue on a physical ground that the sub- 
correction term to the area law for gapped systems typi- 
cally scales as |Ap~^.[47] For such systems, the maximal 
number of topologically protected qubits is bounded by 

- Bounds for more generic systems 

In general, one cannot expect the leading terms of the 
entanglement entropy to be canceled out by choosing an 
appropriate set of subsystems, especially for critical sys- 
tems and ground states of nonlocal Hamiltonian. We 
show that, even for such generic systems, a nontrivial 
tradeoff bound can be obtained. Using the subadditivity 
of entropy, i.e. I{A : B) = S{A) + S{B) - S{AB) > 0, 
Eq|6] yields the following inequality: 



k<J2six,), \X,\<d 



(8) 



for a quantum code with a code distance d and a number 
of encoded qubits k — logA^, where {Xi} is a partition 
of the system. Dividing both sides of the inequality by 
the number of particles n, we find that the rate of a 
quantum error correcting code - is bounded by the aver- 
age entanglement entropy per volume over any partitions 
{Xi}, \Xi\ < d. Our result shows that studying the en- 
tanglement properties of a quantum error correcting code 
is a relevant problem for understanding its fundamental 
limit. In particular, Eq[8] gives a necessary condition for 
a quantum error correcting code to have a nonvanishing 
rate: its average entanglement entropy over subsystems 
smaller than the code distance must satisfy a strict vol- 
ume law. 

EqJ8] can be used to obtain a tradeoff bound for quan- 
tum codes satisfying a subvolume law of entanglement 
entropy. If the entanglement entropy satisfies a subvol- 
ume law, 

S{A) = 0{\An, 

we can obtain the following tradeoff bound; 



kd^-°' = 0{n). 



(9) 



There are at least two important differences between Eq 
and Eq|9] On one hand, Eqj9] is more general than Eq 1 
in that it does not require any structure about the parent 
Hamiltonian. On the other hand, Eqj9] provides a weaker 
tradeoff bound than EqjTJdoes. 

- Outlook 

MED provides a powerful numerical framework 
for studying finite-temperature quantum many-body 
systems. [131 HI] By extending the observation of Hast- 
ings and Poulin, we have shown that MED can be also 
used to obtain a fundamental limit on the information 
storage capacity of a quantum many-body system. In 
particular, we have shown that there is a fundamental 
reason behind the existence of the topological entangle- 
ment entropy. After all, a state with a sufficiently small 



topological entanglement entropy cannot have another 
state that is locally indistinguishable from the original 
state. It is quite remarkable that such a tight tradeoff 
bound can be obtained simply from the property of the 
wavefunction alone. There are several interesting physi- 
cal systems for which our approach might be useful. 

For example, our approach should be easily gener- 
alizable to finite temperature systems. In such sys- 
tems, the topological entanglement entropy is interpreted 
as an order parameter characterizing the phase. [5^55] 
While the technical details must be worked out, we 
can have an intuitive explanation on why those terms 
must be ^2(1). Some of these models become a self- 
correcting classical/quantum memory depending on the 
temperature. [551 [57] Since the "logical operators" that 
are associated to these systems are highly nonlocal, two 
distinct states cannot be distinguished via any local op- 
eration. If the topological entanglement entropy becomes 
sufficiently small, the number of locally indistinguishable 
states is bounded by a small number, which would give 
rise to a contradiction. It is well-known that EqJ4| is 
modified if the subsystem contains a quasiparticle with 
a nontrivial topological charge. [5] It would be interesting 
to extend our analysis to such settings as well. 

We conclude with a remark that our result can be 
used as a rigorous tool for proving a distinctiveness of 
different quantum many-body phases. Given a quan- 
tum many-body system on a torus with N topologi- 
cally protected ground states, it cannot be adiabati- 
cally connected to any state with a topological entan- 
glement entropy strictly lower than °| . This is due to 
the fact that topologically protected states remain to be 
so under an adiabatic evolution. [37l [58] Therefore, the 
lower bound for the topological entanglement entropy, 
i.e. Eqlsl remains stable with a small correction that 
vanishes in the thermodynamic limit. 
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Derivation of EqJ5] under a power law correction 

In the main text of this Letter, we have derived an 
inequality between the topological entanglement entropy 
7 and the ground state degeneracy iV on a torus. There 
are several assumptions that we have implicitly assumed 
in the analysis. First, we have assumed that the entan- 
glement entropy of a region that can be contained in a 
radius of cL for a sufficiently large value of < c < 1 is 
expressed as Eq|4] However, we did not specify the value 
of c explicitly. We have also assumed that the corrections 
to the area law and the local indistinguishability prop- 
erty decays exponentially in the system size. Here we 
show that these conditions can be relaxed significantly 
in general. 

More precisely, we obtain an inequality analogous to 
Eq[5] under the following set of assumptions. 

• The system satisfies a (cL, e)-TQO condition with 
< c < 1, with e = 0{j^), ai > 0. 

• The number of locally indistinguishable states sat- 
isfying the TQO condition is N . 



• The correction to the area law decays algebraically, 
I.e. S{A) = a\dA\ - 7 + O(p^), ^2 > 0. 

• The Hilbert space dimension of the local particles 
is d. 

We begin by bounding a difference between the en- 
tanglement entropy of a single wavefunction, say \^i)^ 
and the entanglement entropy of a mixed state p = 
F Si=i IV'i) (V'il- Without loss of generality, we shall 
assume that Eq|4] is the entanglement entropy formula 
for IV'i)- Denoting the entanglement entropy of the sin- 
gle wavefunction as S^{A) and the entanglement entropy 
of the mixed state p as Sp{A), we obtain the following 
bound. 

Lemma 1. Let A he a region that can he contained in a 
hall of radius cL. 



Also, 



\S\A) ^SP{A)\<e{d\A\- log e) 



(10) 



This result follows from a simple application of the 
Fannes inequality, see EqlT] 

As we did in the main text, consider a sequence of 
subsystems Ai, Bi, and Ci such that AiBiCi — Ai+iBi+i. 
By choosing the subsystems AiBi and BiCi to be small 
enough to be contained in a ball of radius cL^ 



logN <SP{AiBi) 



i=l 



SP(B,C,) - SP[B,) 



(11) 



< S\A^B^) + Y, S\B,a) ~ S'iB,) + eLA, 

(12) 



i=l 



where eLA is the approximation error from Eq |10| One 
can choose the size of the subsystems AiBi and BiCi to 
be Q{c^L^), in which case n becomes n = ©(-^)- There- 
fore, 



eLA = 0{edL^). 



(13) 



S\A,B,) + J2 S\B,Ci) - S\B,) = 27 + 0{—^ 



2(1 + Q2) 






L2"2 



(14) 



Combining these bounds together, for a constant c and 
d, we conclude that 



log iV < 27 + 0{L'-"^ ) + 0{L 



(15) 



Therefore, for a sufficiently large system size, our main 
result holds for any ax > 2, a2 > 0. The signifi- 
cance of this bound comes from the fact that a conser- 
vative estimate for the area law correction term gives 
02 — ^ • HZ] The preceding analysis shows the robust- 
ness of our bound even under such a conservative error 
estimate. Even under an algebraic correction to the ideal 
wavefunction, the inequality remains intact. 

We would also like to explain the intuition behind why 
a factor of 2 appears in front of 7. It should be clear from 
the construction that the area terms cancel out. There- 
fore, one only needs to be concerned about the topology 
of the subsystems AiBi, BiCi, and Bi. The general idea 
is that for the construction in FIGJl] AiBi and BiCi are 
simply connected, yet Bi is a union of two simply con- 
nected subsystems. Hence the topological contribution 
must be 7 for AiBi and i?iCi, and 27 for Bi. Similarly, 
B2C2 gives a contribution of 7, while B2 gives 27. B3C3 
is simply connected, so it gives a contribution of 7. B^ 
is an annulus, hence it gives 27. Combining these results 
together, the topological contribution to our bound be- 
comes 27. A similar analysis can be carried out for other 
choice of subsystems as well. What is important is that 
Eq[TT] holds for any choice of subsystems. Therefore, one 
can always optimize over different set of subsystems to 
obtain the best bound. 



